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ABSTRACT 

The origin of the irregular satellites of the giant planets has been long debated since 
their discovery. Their dynamical features argue against an in-situ formation suggesting 
they are captured bodies, yet there is no global consensus on the physical process at 
the basis of their capture. In this paper we explore the collisional capture scenario, 
where the actual satellites originated from impacts occurred within Saturn's influence 
sphere. By modeling the inverse capture problem, we estimated the families of orbits 
of the possible parent bodies and the specific impulse needed for their capture. The 
orbits of these putative parent bodies are compared to those of the minor bodies of the 
outer Solar System to outline their possible region of formation. Finally we tested the 
collisional capture hypothesis on Phoebe by taking advantage of the data supplied by 
Cassini on its major crater, Jason. Our results presented a realistic range of solutions 
matching the observational and dynamical data. 

Key words: planets and satellites: formation, Saturn, irregular satellites - methods: 
numerical, N-Body simulations - celestial mechanics 



1 INTRODUCTION 

In our previous paper (Turrini et al., accepted, hereafter 
Paper I) on the subject of Saturn's irregular satellites we 
explored the dynamical and collisional structure of the sys- 
tem to shed new light on its past history. The results we 
obtained strongly suggested that the system of irregular 
satellites of Saturn went subject to a strong dynamical, and 
probably also collisional, sculpting during its lifetime. The 
actual structure of the system appeared to be very likely 
less representative than previously thought of the original, 
post-capture one. Nevertheless, through our results we have 
been able to draw a coherent picture of the possible past 
evolution of the system. 

Before proceeding further, we would like to briefly sum- 
marise the main results of our previous work: 

• we found confirmation for the resonant motion in the 
Kozai regime for Ijiraq and, partially, Kiviuq; 

• our numerical simulations showed strong evidences of 
the presence of chaos in the system, with about two thirds 
of the satellites being affected to some extent; 

• we verified that the influence of the Great Inequality 
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between Saturn and Jupiter plays a central role in shaping 
the dynamical behaviours of the irregular satellites; 

• we found that the influence of the two outermost major 
satellites (Titan and Iapetus) alters the evolution of the sys- 
tem, having stabilising effects on some satellites and desta- 
bilising ones on others; 

• we verified that the present orbital structure is long- 
lived against collisions but there are indications that in the 
past a more intense collisional activity could have taken 
place; 

• we interpreted the absence of prograde irregular satel- 
lites in the region comprised between 11.22 x 10 6 km and 
14.96 x 10 6 km as a by-product of the sweeping effect of 
Phoebe; 

• similarly, we suggest that the absence of retrograde 
satellites in the same radial region could be due to the same 
mechanism, even if we could not rule out the possibility of 
this being a primordial heritage; 

• by applying the HCM algorithm (IZappala et al. 111990 



1994 1 we found two candidate families between the prograde 



and six between the retrograde satellites, with some of the 
families merging into bigger groups we called clusters. 

The framework of our first article was our project to study 
the origin of the irregular satellites, using Saturn system 
as our case study to take advantage of the unprecedented 
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data supplied by the Cassini mission during its only flyby of 
Phoebe. In particular, our interest concentrated on investi- 
gating the viability and the consequences of the collisional 
capture process, originally proposed by |Colombo fc Franklin 
( 1971 1, which could supply a natural way to explain the for- 



mation of the systems of irregular satellites of both gas and 
ice giant planets. With the deeper insight of the dynamical 
and collisional nature of the system of irregular satellites 
of Saturn obtained through our previous results, we investi- 
gated the collisional capture process with the aims of con- 
straining the characteristics of the impacts and of searching 
for information on the possible pre-capture orbits. 
The work to be presented is organised as follows: 

• in section [2] we will describe the collisional capture sce- 
nario and the dynamical strategy we devised to explore its 
implications; 

• in section [3] we will compare the pre-capture orbits ob- 
tained through the numerical code we developed with the 
ones of the minor bodies in the outer Solar System; 

• in section [4] we will search for evidences supporting 
the collisional capture hypothesis in the data supplied by 
Cassini; 

• finally, in section [5] we will put together all our results 
to draw a global picture of the origin of Saturn's system of 
irregular satellites. 



2 THE COLLISIONAL CAPTURE SCENARIO 

Our interest in studying the capture mechanisms of ir- 
regular satellites concentrated on the collisional scenario 



originally proposed by Colombo & Franklin (19711, which 



seems to be the only one capable to explain in a natural 
way the existence of the systems of irregular satellites of 
both gas and ice giants. 

The ideal, ground-based strategy to investigate this capture 
scenario would be to identify the possible collisional families 
existing in the system through the comparative study of 
their mean orbital elements and their surface compositions. 
In principle, the masses, the orbits and the composition of 
the parent bodies, once estimated, could make it possible 
to identify their nature and their region of formation. Even 
more ideal it would be to search for features proving the 
occurrence of relevant or disruptive impacts on the surface 
of the satellites through space-based observations from 
probes performing close fly-bys. 

The results we presented in our previous paper showed that 
these ideal conditions are difficult to meet in the case of 
the Saturn system. The collisional history of its irregular 
satellites appears to have been eventful and the actual 
structure of the system is probably not representative of 
the original, post-capture one. Moreover, the combined 
effects of secondary collisions and of the perturbed, chaotic 
dynamical evolution of some of the satellites could have 
contributed in mixing and altering the possible families 
generated by catastrophically disruptive events, thus in- 
creasing the difficulty of the task of identifying them. 
In addition to these problems, the distance and the small 
size of the irregular satellites make it difficult to observe 
them and obtain detailed spectral information with Earth- 
based telescopes. In general, the only observational data 
that could be gathered are the information-poor colour 



indexes. Even for the biggest and most luminous irregular 
satellites, the spectral range achievable is quite limited and 
it is therefore difficult to trace back their compositions. 
In Saturn's case, however, one of the ideal conditions for 
the study of the collisional capture scenario can be met. 
Thanks to the Cassini spacecraft, which performed the 
first and to now only observations of an irregular satellite 
during a close fly-by, detailed information on the physical 
and compositional nature of Phoebe are available. This is 
the reason behind our choice of the Saturn system as our 
case study. 



2.1 The inverse capture problem: dynamical 
model and numerical code 

Since realistic simulations of the collisional capture process 
are to now an unreachable goal due to the limited data avail- 
able, the arbitrary initial conditions and the degeneracy of 
the solutions, we devised the following strategy to investi- 
gate the problem. Instead of exploring directly the collisional 
capture, we concentrated on the inverse problem (i.e. the 
orbital unbinding from the host planet) taking advantage 
of the time reversibility of the equations of motion. The 
dynamical model we devised employs the specific impulse 
needed to obtain the passage from a bound planetocentric 
orbit to a bound heliocentric orbit as its only parameter and 
is designed to search for the possible related solutions. Our 
choice to employ the specific impulse as the free parameter 
aimed at minimising the dependence of the solutions from 
the characteristics of the involved bodies and to obtain in- 
formation as general as possible. The assumption underlying 
the model is that the capture of an irregular satellite should 
have been due to a single, predominant impact event which 
we could model through the application of a delta-shaped 
specific impulse (i.e. an instantaneous change in velocity). 
We numerically implemented the model in the code MSSCC 
(Modelling Software for Satellite Collisional Capture), based 
on the patched conies approximation, which solving algo- 
rithm addresses the inverse problem in the following way: 

• given the initial planetocentric orbital elements of the 
satellite, its (f,v) planetocentric vectors are computed; 

• a spherical sampling of vectors dvl is created using the 
fixed modulus dv and the angular sampling steps A9 and 
Aip specified in input by the user; 

• the specific impulse vectors dvl are summed to the v 
vector and the new pairs of conjugate variables (fi,Vi ) are 
computed with fl = r ; 

• converting the (fl,v~l ) pairs to orbital elements, the or- 
bits which are hyperbolic in the planetocentric frame are 
kept while the others are discarded; 

• the selected orbits are reverted to the heliocentric frame 
and only the elliptic, prograde ones are saved as candidate 
primordial orbits for the parent bodies. 

The condition on the inclination of the candidate orbits (i.e. 
they must be prograde) is imposed to assure the physical 
correctness of the solutions found, since the formation of 
retrograde bodies from the Solar Nebula is not plausible. 
If specified by the user, the code can loop the algorithm while 
varying the orbital position of both the satellite around the 
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planet and the planet around the Sun to avoid selection ef- 
fects due to the choice of particular initial conditions. The 
code also implements the possibility to filter the candidate 
orbits over a selected e m ax value of the eccentricity, in order 
to control the range of physical conditions under which the 
capture event took place. 

An open issue with this approach is the physical meaning of 
the specific impulse dv used by the algorithm. To the first 
order and for ideally, fully inelastic collisions the interpre- 
tation is straightforward, since it represents the change in 
velocity due to the collision of the parent body with the 
projectile. The real case is however complicated by the fact 
that impacts always involve some level of energy dissipa- 
tion due to the generation of heat, the formation of craters 
and fracture lines in the target body and the excavation of 
fragments. Moreover, things are worse in case of catastroph- 
ically disruptive events. While for the ejection of fragments 
the total excavated mass is usually small in comparison to 
the mass of the target and the contribution to the change 
in momentum is limited, in case of catastrophic disruption 
the mass of the collisional shards could be relevant and the 
ejection speeds, estimated from laboratory experiments and 
numerical simulations with hydrocodes to be of the order of 
100 m/s (Benz & Aspaugh 19991, could represent a signifi- 
cant fraction of the specific impulse imposed by the collision. 
While this problem in general has to be dealt with on a case 
by case basis, we will treat it in major detail for Phoebe in 
section [4] 



2.2 Settings and results of the dynamical analysis 

The configuration we used for our investigation of the inverse 
capture process of Saturn's irregular satellites has been: 

• sampling step for the dvl vectors in the x — y plane: 5° 

• sampling step for the dvl vectors in the z direction: 5° 

• sampling step for the satellite mean anomaly: 5° 

• sampling step for the planetary mean anomaly: 5° 

• limiting value for heliocentric eccentricity: 0.9 

Since our model employed the specific impulse as its free pa- 
rameter, we were not constrained by the size of the bodies in 
the choice of the ones to be used as our reference cases. We 
chose four irregular satellites (two prograde and two retro- 
grade) representative of the main dynamical configurations 
for our investigation: Albiorix, Siarnaq, Phoebe and Mundil- 
fari. As their orbits, we used the mean orbital elements we 
presented in Paper I. Phoebe was a natural choice as one of 
our case studies since our simulations showed that the sec- 
ular variations of its orbit are quite regular and limited (see 
Paper I), implying a better preservation of its primordial 
dynamics. The results we obtained are presented from fig. 
[I] to fig. |16| where we showed, for increasing values of the 
specific impulse, the evolution of the solutions to the inverse 
problem in the: 
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respectively 10 4 , 3 x 10 4 , 10 5 , 2 x 10 5 , 6 x 10 s , 1.2 x 10 6 solu- 
tions. These values have been chosen to represent the main 
features of the evolution of the solutions and range from the 
threshold values for which the first solutions appeared to 
the values for which the different families of solutions over- 
lapped into a single continuum. We found that the minimum 
change in velocity to be applied in order to capture the ir- 
regular satellites lay in general between 450 m/s and 500 
m/s, with Phoebe being a separate case requiring a value 
(« 650 m/s) about 30% higher. 

While the evolution of the solutions of each satellite showed 
its own peculiar features, there were some general conclu- 
sions we could draw. First, we noticed that the solutions 
were grouped into distinct and initially well separated fami- 
lies (see top panels of figs. Tfl6 I which, for increasing values 
of the specific impulse applied, tended to overlap to finally 
merge into a more or less continuous distribution (see bot- 
tom panels of figs. 1)16 1. Second, these families of solutions 
were related to preferential direction of approach, as is par- 
ticularly evident from the top panels of figs, [I] [5] [9] and |13| 
This occurrence is strictly related to the existence of prefer- 
ential directions of impulse change (i.e. impact geometry) as 
is apparent from the plots of figs. |4|[8| |12| and |16| Third, the 
first solutions to appear for the lowest values of specific im- 
pulse applied were characterised by high eccentricity values 
(generally 0.3 < e < 0.7, except in Mundilfari's case where 
0.2 < e < 0.5) and semimajor axis spread between 5 — 20 
AU (5 — 30 AU in Albiorix's case), a behaviour contrary to 
what was generally believed to occur (e.g. the reservoir of 
parent bodies of the irregular satellites being located near 
their host planets). 

As a rule of the thumb, prograde and retrograde satellites 
formed two distinct groups in terms of the solutions to the 
inverse capture problem. If we compare the cases with the 
same number of solutions, the prograde cases cover a wider 
region of phase space than retrograde cases. The range of 
inclination values covered by the former is at least one 
third higher than the one of the latter (0° < i < 12° vs. 
0° < i < 8°, see figs. [5J [6] [ToJ [U] and figs.[3j0[nj[T5]), the 
semimajor axis of the solutions reaches regions more distant 
from Saturn (see figs.[Tj[5"||9| |13[ ) and the eccentricity reaches 
both higher and lower values (e < 0.2 and e > 0.7 respec- 
tively, see figs. [5] §[13]). 

The solutions related to the orbital region inside Saturn's 
orbit are the most widespread in inclination and acquire the 
highest inclination values (see bottom panels of figs. [2] [6] 



101 \14\ . The solutions of Albiorix, Siarnaq and Phoebe ini- 
tially avoid low inclination values (i < 3°) for those orbits 
coming from the outer part of the Solar System with re- 
spect to Saturn's orbit (see top panels of figs. |2l 
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inclination values are acquired for higher values of the spe- 
cific impulse applied. Mundilfari had an opposite behaviour 



The values of the specific impulse employed for each satellite 
(reported in the captions of the associated figures) produced 



(see top panels of fig. 14 1 and avoided low inclination values 
for those orbits with semimajor axis located inside Saturn's 
orbit. There was no apparent correlation between the ec- 
centricity and the inclination values of the solutions, with 
two major exceptions for those solutions associated respec- 
tively with the lowest (e < 0.2) and the highest (generally 
e > 0.5) eccentricity values. The former case is related to 
the outermost solutions, which avoid low inclination values 
(i < 2° — 3°) and tend to cluster between 3° < i < 5° (see 
bottom panels of figs. |3j [7| |11| |15| ), Albiorix had a second 
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Figure 1. Evolution of the solutions to the inverse capture problem obtained for Albiorix for increasing values of the specific impulse 
applied and projected in the heliocentric a — e plane. From top left to bottom right, the values assumed by the specific impulse are 
445, 460, 520, 580, 760, 880 m/s. Distances are expressed in AU. 



tail in the range 5° < i < 8° (see fig.[3| bottom panels). The 
latter case is instead related to those orbits spatially located 
near Saturn and appears in the figures associated to the pro- 
grade satellites Albiorix and Siarnaq (see fig. [3] and fig. [7J 
bottom panels). For those solutions we observed an inverse 
correlation between the inclination and the eccentricity (i.e. 
the lower the eccentricity, the higher the inclination). 
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Figure 2. Evolution of the solutions to the inverse capture problem obtained for Albiorix for increasing values of the specific impulse 
applied and projected in the heliocentric a — i plane. From top left to bottom right, the values assumed by the specific impulse are 
445, 460, 520, 580, 760, 880 m/s. Distances are expressed in AU while angles are expressed in degrees. 
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Figure 3. Evolution of the solutions to the inverse capture problem obtained for Albiorix for increasing values of the specific impulse 
applied and projected in the heliocentric e — i plane. From top left to bottom right, the values assumed by the specific impulse are 
445, 460, 520, 580, 760, 880 m/s. Angles are expressed in degrees. 
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Figure 4. Specific impulse applied to Albiorix to solve the inverse capture problem projected in the v x — v y plane. From top left to 
bottom right, the values assumed by the specific impulse are 445, 460, 520, 580, 760, 880 m/s. Velocities in the graphs are expressed in 
km/s. 
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Figure 5. Evolution of the solutions to the inverse capture problem obtained for Siarnaq for increasing values of the specific impulse 
applied and projected in the heliocentric a — e plane. From top left to bottom right, the values assumed by the specific impulse are 
490, 500, 540, 580, 700, 800 m/s. Distances are expressed in AU. 



Saturn Irregular Satellites - Origin 9 




Figure 6. Evolution of the solutions to the inverse capture problem obtained for Siarnaq for increasing values of the specific impulse 
applied and projected in the heliocentric a — i plane. From top left to bottom right, the values assumed by the specific impulse are 
490, 500, 540, 580, 700, 800 m/s. Distances are expressed in AU while angles are expressed in degrees. 
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Figure 7. Evolution of the solutions to the inverse capture problem obtained for Siarnaq for increasing values of the specific impulse 
applied and projected in the heliocentric e — i plane. From top left to bottom right, the values assumed by the specific impulse are 
490, 500, 540, 580, 700, 800 m/s. Angles are expressed in degrees. 
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Figure 8. Specific impulse applied to Siarnaq to solve the inverse capture problem projected in the v x — v y plane. From top left to 
bottom right, the values assumed by the specific impulse are 490, 500, 540, 580, 700, 800 m/s. Velocities in the graphs are expressed in 
km/s. 
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Figure 9. Evolution of the solutions to the inverse capture problem obtained for Phoebe for increasing values of the specific impulse 
applied and projected in the heliocentric a — e plane. From top left to bottom right, the values assumed by the specific impulse are 
640, 655, 685, 720, 800, 880 m/s. Distances are expressed in AU. 



Saturn Irregular Satellites - Origin 13 




Figure 10. Evolution of the solutions to the inverse capture problem obtained for Phoebe for increasing values of the specific impulse 
applied and projected in the heliocentric a — i plane. From top left to bottom right, the values assumed by the specific impulse are 
640, 655, 685, 720, 800, 880 m/s. Distances are expressed in AU while angles are expressed in degrees. 
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Figure 11. Evolution of the solutions to the inverse capture problem obtained for Phoebe for increasing values of the specific impulse 
applied and projected in the heliocentric e — i plane. From top left to bottom right, the values assumed by the specific impulse are 
640, 655, 685, 720, 800, 880 m/s. Angles are expressed in degrees. 
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Figure 12. Specific impulse applied to Phoebe to solve the inverse capture problem projected in the v x — v y plane. From top left to 
bottom right, the values assumed by the specific impulse are 640, 655, 685, 720, 800, 880 m/s. Velocities in the graphs are expressed in 
km/s. 
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Figure 13. Evolution of the solutions to the inverse capture problem obtained for Mundilfari for increasing values of the specific impulse 
applied and projected in the heliocentric a — e plane. From top left to bottom right, the values assumed by the specific impulse are 
515, 530, 560, 590, 670, 750 m/s. Distances are expressed in AU. 
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Figure 14. Evolution of the solutions to the inverse capture problem obtained for Mundilfari for increasing values of the specific impulse 
applied and projected in the heliocentric a — i plane. From top left to bottom right, the values assumed by the specific impulse are 
515, 530, 560, 590, 670, 750 m/s. Distances are expressed in AU while angles are expressed in degrees. 
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Figure 15. Evolution of the solutions to the inverse capture problem obtained for Mundilfari for increasing values of the specific impulse 
applied and projected in the heliocentric e — i plane. From top left to bottom right, the values assumed by the specific impulse are 
515, 530, 560, 590, 670, 750 m/s. Angles are expressed in degrees. 
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3 PARENT BODIES OF THE IRREGULAR 
SATELLITES 

Once we completed the first part of our investigation of 
the collisional capture problem, we proceeded to the second 
stage: the search for the primordial populations the parent 
bodies of Saturn's irregular satellites could have originated 
from. To reach such an ambitious goal, we compared the 
dynamical solutions we found and the observational data 
about Saturn's irregular satellites to all the available dy- 
namical and observational data on the populations of minor 
bodies in the outer Solar System. Moreover, to test the ca- 
pability of the collisional mechanism to explain the existence 
of the irregular satellites in the scenario outlined by the Nice 
model, we extended the comparison to the bodies crossing 
the orbit of Saturn in the original simulations described in 
Gomes et aT] | |2005| |, [Morbidelli et aI7] ( |2005p and |Tsiganis~et| 
al. (20051. The dynamical data on the minor bodies of the 



outer Solar System have been obtained by the catalogues 
of the IAU Minor Planet CenteiQand referred to the orbits 
determined at November 2006. 

In figure |17[ |18| and [19] we show the solutions we found for 
the four Saturn's irregular satellites we used as our case 
study (i.e. Albiorix and Siarnaq for the prograde group 
and Mundilfari and Phoebe for the retrograde one). Super- 
imposed to the solutions we plotted the families of minor 
bodies of the outer Solar System: comets, trans-neptunian 
objects and Scattered Disk objects (respectively TNOs and 
SDOs) and Centaurs. We also plotted the synthetic family 



of Saturn Crossers found by Gomes et al. (20051, Morbidelli 



et al. ( 2005 I and Tsiganis et al. ( 2005 I in their studies of the 



formation of the orbital configuration of the Solar System 
following a dynamical rearrangement leading to the Late 
Heavy Bombardment 

As can be seen from the figures, three of the four fami- 
lies of minor bodies considered (Centaurs, comets and Sat- 
urn Crossers, TNOs and SDOs being dynamically decoupled 
from Saturn) overlap to various extents the region of phase 
space occupied by the solutions to the inverse capture prob- 
lem: 

• comets match the solutions for eccentricities in the 
range 0.1 < e < 0.6 and in the radial region 6AU < a < 
10 AU ; the retrograde cases require i < 8° , while for pro- 
grade cases inclination values i < 12° are required; 

• Centaurs partially overlap the solutions for semimajor 
axis ac and eccentricity ec values satisfying the relationship 
ac(l — ec) < a-s with as being Saturn's semimajor axis. The 
constraints on inclination are the same as found for comets; 

• Saturn Crossers partially overlap the space occupied by 
the solutions in the radial interval 7AU < a < 70AU for 
orbits having e > 0.1; again, the constraints on inclination 
are the same as found for comets. 

Once we completed the dynamical comparison of the solu- 
tions to the minor bodies of interest, we proceeded in verify- 
ing if the match we found extended also to the spectropho- 
tometric data (i.e. the colour indexes). 



Various authors (Grav et al. 2003 Grav & Holman 2004 



Buratti et al.|2005 I already performed ground-based photo- 
metric observations of the irregular satellites of Saturn and 
compared their colours with the ones of Centaurs, TNOs and 



2006 



recently with cometary nuclei and dead comets (Sheppard 
Jewitt|[2008 1. The results of such comparisons can be 



summarised as follows: 

• irregular satellites present colours comparable with 
those of grey TNOs and of Centaurs; 

• irregular satellites lack the so-called ultrared matter 
characteristic of the spectra of many TNOs and Centaurs; 

• irregular satellites have colours distributions similar to 
those of cometary nuclei. 

Taking advantage of the results of our dynamical compar- 
ison, we devised the following strategy: instead of compar- 
ing the whole populations of Centaurs, comets and irregu- 
lar satellites as previously done, we decided to compare the 
observational data of the irregular satellites with those of 
the subsets of Centaurs and comets matching the dynami- 
cal range of the solutions to the inverse problem. The aim 
behind our strategy was to perform a targeted comparison to 
minimise the effects of extraneous or of dynamically incom- 
patible bodies, which could invalidate the global picture. 
We proceeded by first gathering all the information avail- 
able in the literature through the major and most updated 
catalogues on the minor bodies in the outer Solar System: 

• MBOSS (Minor Bodies of the Outer Solar SysterrQ 

• 2MASS Comet Surve/l 

• COMETS II catalogue}^ 

• 2MS (Meudon Multicolour SurvejQ) 

The collected data were then cross-correlated with the dy- 
namical information from the IAU Minor Planet Center's 
catalogues, to link each body to its orbits. 
The observational data on Saturn's irregular satellites were 



taken from Grav et al. ( 2003 1 (BVRI bands) and from Grav 
|fc Holman| ( |2004[ ) (JHK bands). 
In considering the observational data we concentrated on 
the colour indexes instead of the spectral data because the 
spectra of irregular satellites and minor bodies of the outer 
Solar System obtained from Earth referred to too small a 
number of bodies and to too small a spectral range to allow 
for significant comparisons. Moreover, we could not take ad- 
vantage of the high-resolution spectra produced by Cassini 
for Phoebe due to the fact that their very resolution and 
their wider spectral range made them unsuitable for a com- 
parison with the range-limited, ground-based data. 



As can be seen from &g.^20\(BVRI bands) and fig. [21] (JHK 
bands, only irregular satellites and comets from 2MASS are 
plotted), by considering la error intervals we found a good 
degree of overlapping between the two groups of bodies in 
the B — V vs V — R plane and in the JHK colours. We found 
a less satisfying match in the planes B — V vs V — I and 
V — R vs V — I due to the V — I colour: however, the V — I 
colours of irregular satellites and minor bodies are compat- 
ible at the 3<r level. 



1 http: //cf a- www. harvard, edu/ iau/ 



Morbidelli, personal communication 



www . sc . eso 


. org/~ohainaut/MB0SS/ 




http : / /www 


. ipac . caltech . edu/2mass/ 


Lamy et al. ( 2004b 
Delsanti et aqi|2Dl)6| 
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Figure 17. Comparison in the a — e plane between solutions to the inverse capture problem (yellow filled diamonds), comets (grey filled 
upper triangles), Centaurs and Scattered Disk objects (violet lower empty triangles), trans-neptunian objects (cyan filled squares) and 
Saturn Crossers (green empty squares). The irregular satellites considered were the prograde satellites Albiorix and Siarnaq and the 
retrograde satellites Mundilfari and Phoebe (clockwise from upper left). The specific impulse values employed for the numerical solutions 
of each irregular satellites are the highest ones we considered, producing 1.2 X 10 6 solutions. Distances are expressed in AU. 



Such results are encouraging since they seem to validate our 
dynamical investigation of the inverse capture problem, but 
we should consider them with care: even if such a matching 
does exist, it exists between the present day populations of 
these families. While the colour indexes of the irregular satel- 
lites should have changed more or less homogeneously due to 
their evolving under constant environmental conditions, we 
cannot say the same for comets and Centaurs. Their dynam- 
ical lifetimes of the order of 10 s years mean that what we 
see today is a cosmogonically young population that could 
be not representative (at least from a compositional point 
of view) of the past ones. Moreover, the different orbital re- 
gions crossed during their lifetimes imply that they evolved 
under different environmental conditions. 



4 COMPARISON WITH CASSINI DATA ON 
PHOEBE 

As we pointed out in the previous section, a direct compari- 
son of the compositional data supplied by the Cassini VIMS 
instrument with the ones available for the other populations 
of minor bodies is hardly possible at the present moment 
because of observational issues. Nevertheless, we found that 
the ISS images of Phoebe could supply us with a number 



of informations extremely precious in the context of our re- 
search: the models of Phoebe's craters produced by the ISS 
images by Giese et al. ( 2006 1 gave us in fact the first physi- 



cal evidence in support of the collisional capture scenario. 
Without information on the dating of the cratering events, 
the interpretation of Phoebe's collisional history is not 
straightforward. We can divide the collisional capture 
scenarios into two main groups: catastrophic and non- 
catastrophic scenarios. In catastrophic scenarios, Phoebe 
would be the largest remnant of a bigger parent body shat- 
tered during the capture event. The craters on the surface 
of the satellite would be the result of the re-accretion of the 
smallest collisional shards produced in the impact. The re- 
sults we presented in Paper I indicate that Phoebe would 
be highly efficient in removing near-by objects in the re- 
gion of phase space the collisional shards would populate, 
explaining both the intense cratering on the surface of the 
satellite and the absence of other remnants of the parent 
body. Non-catastrophic scenarios would imply that Phoebe 
and its parent body share a similar mass, modified only by 
the fragmentation and excavation effects of post-capture im- 
pacts. Both kinds of scenario suggest that the dynamical 
evolution of Phoebe could have suffered major changes in 
the past as a consequence of impact events. 
Catastrophic scenarios are difficult to model because of the 
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Figure 18. Comparison in the a — i plane between solutions to the inverse capture problem (yellow filled diamonds), comets (grey filled 
upper triangles), Centaurs and Scattered Disk objects (violet lower empty triangles), trans-neptunian objects (cyan filled squares) and 
Saturn Crossers (green empty squares). The irregular satellites considered were the prograde satellites Albiorix and Siarnaq and the 
retrograde satellites Mundilfari and Phoebe (clockwise from upper left). The specific impulse values employed for the numerical solutions 
of each irregular satellites are the highest ones we considered, producing 1.2 X 10 e solutions. Distances are expressed in AU, angles are 
expressed in degrees. 



higher degeneracy in the initial conditions (e.g. the target 
and the impactor mass values) than non-catastrophic ones. 
We therefore decided to investigate the non-catastrophic 
scenarios and verify if the possible impacts which produced 
the major craters, under the right conditions, could have 
applied a specific impulse of the order of the ones previ- 
ously estimated and thus explain the capture of the satel- 
lite. In principle, once in possession of information on the 
morphology, the depth and the diameter of a crater, it is pos- 
sible to estimate the characteristics of the impact and the 
projectile by inverting the empirical formula bylSchmidt & 



Housen (19871. Following Nesvorny et al. (2003 



and 



Pier- 



azzo fc Melosh| ( |2000| ) , we modified the formula to account 



for the dependence on the vertical component of impact ve- 



locity ViCosa where a is the impact zenith angle (Gault & 
|Wedekind|1978t : 



V = 0.13 



Pt 



Pt 



(t 



(1) 



where all quantities have to be evaluated in cgs units and: 

• rrii, pi are the mass and density of the impactor 

• p t is the density of the target 

• Vi is the velocity of the impactor 



• a is the incidence angle measured from the zenith 

• g is the surface gravity of the target 

• V is the volume of the crater 

In general, however, due to the degeneracy of the initial con- 
ditions it's not possible to obtain a deterministic solution to 
the problem. Moreover, Phoebe's case is complicated by the 
fact that Cassini images revealed the existence of at least 7 
major craters, distributed inhomogeneously on the surface 
of the satellite and with diameters greater than 50 km (in 
two cases, Jason and Eurytus craters, the diameter is about 
100 km). We concentrated our attention on Jason crater, for 
which the resolution and the completeness of Cassini data 



allowed the development of a structure model (Giese et al. 
2006[). Th e analysis of Hylas crater (30 km diameter, |Giese et 
al. ( 2006 1) revealed a depth to diameter ratio 1 : 6, while for 
the minor craters this value is about 1 : 5. Since equation[l]is 
strictly appliable only to geometrically simple, bowl-shaped 
craters while Jason crater shows a complex structure, made 
more difficult to interpret by the presence of landslides and 
of a secondary crater superimposed, we relied on the follow- 
ing approximation. We assumed a bowl-shaped structure 
using the diameter estimated by Cassini measures (101 km) 
and a depth to diameter ratio of 1 : 7. This value was chosen 
to make a conservative estimate of Jason's volume and to 
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Figure 19. Comparison in the e — i plane between solutions to the inverse capture problem (yellow filled diamonds), comets (grey filled 
upper triangles), Centaurs and Scattered Disk objects (violet lower empty triangles), trans-neptunian objects (cyan filled squares) and 
Saturn Crossers (green empty squares). The irregular satellites considered were the prograde satellites Albiorix and Siarnaq and the 
retrograde satellites Mundilfari and Phoebe (clockwise from upper left). The specific impulse values employed for the numerical solutions 
of each irregular satellites are the highest ones we considered, producing 1.2 X 10 6 solutions. Angles are expressed in degrees. 



take into account its complex morphology. The depth esti- 
mated this way (14.43 km) is compatible with the range of 
values (13 — 20 km) reported by Giese et al. (2006). Through 
the formula 



In the analysis we used the value of 5 cm/a 2 for Phoebe's 



V = nh x 



h 2 
6 



(2) 



with a 



we obtained an estimated volume of 5.94 x 10 
total excavated mass equal to 9.68 x 10 19 g if we assume its 
density equal to Phoebe's mean density (1.63 g/cm 3 



Porco 



et al. ( 2005 1). Since Phoebe's mass, estimated by using the 
mean radius (106.6 km, Porco et al. (2005)) and the mean 
density supplied by Cassini measures, is 8.26 x 10 21 g, we 
can easily verify that, for ejection velocities of the order of 



100 m/s as estimated by Bcnz & Aspaugh (19991, the total 



contribution to the change in momentum due to the ex- 
cavated fragments is negligible with respect to the specific 
impulse necessary for Phoebe's capture. 
By inverting eq. [T] we obtain a relationship linking the im- 
pact speed to the characteristics of both impact and crater: 



V 



0.13 



(cos a) 1 



(3) 



surface gravity ( Porco et al.|2005 | and we attributed to the 
density of the impactor the same value as Phoebe's mean 
density. The free parameters in our analysis were the mass 
of the impactor rrii and the incidence angle a. Note that the 
values of both densities and of the surface gravity, due to 
their range of variation, do not change the order of magni- 
tude of the results. 

Before proceeding with the analysis we want to spend a few 
words on the link between the impact velocity and the spe- 
cific impulse necessary to capture the satellite. The impact 
speed can be related to the impulse transmitted to Phoebe 
from the collision by considering the motion of both im- 
pactor and target in a reference frame co-moving with the 
latter before the impact. In this reference frame the target 
is initially at rest while the impactor moves at the relative 
speed Vrei which, save for gravitational focusing effects, is 
the same as the impact speed. The gravitational effects can 
be evaluated by the escape velocities of the two bodies, once 
their masses are known or have been fixed, and we can con- 
sider them through a corrective parameter: for Phoebe, the 
escape speed is of the order of lOOm/s ( Porco et al.|2005 l. 
In general, for real, inelastic collisions we cannot solve the 
equations for the conservation of the linear momentum since 
we don't know the amount of momentum carried away by 
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Figure 20. Colour-colour diagrams in the BVRI spectral region 
of irregular satellites (black symbols) and minor bodies from the 
catalogues considered whose orbits match the dynamical region 
described by the solutions to the inverse capture problem (grey 
symbols). The error bars represent the lcr confidence level of the 
colour indexes, where available. 



Figure 21. Colour— colour diagrams in the JHK spectral region 
of irregular satellites (black symbols) and comets from 2MASS 
catalogue whose orbits match the dynamical region described by 
the solutions to the inverse capture problem (grey symbols). The 
error bars represent the lcr confidence level of the colour indexes, 
where available. 



the impactor (i.e. the amount of momentum transferred) and 
the general efficiency of the process, which depends on the 
physical characteristics of the two bodies. As a consequence, 
to give a rough estimate of the plausibility of the collisional 
capture mechanism, we initially approximated the collision 
as perfectly elastic and ignored the effects of the excava- 
tion and ejection of fragments (which, for the range of mass 
and velocity values here considered and the typical ejection 



speeds, would amount to less than 0.1% of the total momen- 
tum). The elastic approximation implies that we expect the 
impactor to ricochet and bounce from the target with a ve- 
locity higher than the escape velocity: the reason behind this 
choice will become clear in the following. In this scenario, 
by varying the free parameters and taking into account the 
gravitational focusing effects (which, depending on the mass 
values considered, could vary between 100—200 m/s), among 
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the possible solutions we found that we could produce the 
change in velocity of ~ 650 m/s necessary to capture Phoebe 
and excavate Jason crater through highly grazing impacts 
(88° < a < 90°). The impactors had mass values varying 
between one tenth of Phoebe's mass and the mass of the ir- 
regular satellite (i.e. 0.8 — 8 x 10 21 g) and were moving at rel- 
ative speeds between lkm/ s < v re \ < 15 km/s (the higher 
velocities relating to the lowest values of mass). For den- 
sity values matching the mean density of Phoebe, that mass 
range would imply the size of the impactor being comprised 
between about 50 km and 100 km. To improve our estimate 
of the plausibility of the collisional capture mechanism, we 
changed our approximation from perfectly elastic to inelastic 
(i.e. the impactor still bouncing with a velocity higher than 
the escape speed) by including in our equations a coefficient 
of restitution equal to 0.5. We obtained the same results 
previously presented with a more strict requirement on the 
impact angle, which now varies in the range 89° < a < 90°. 
The range of variations of the impact angle is at the basis 
of our assumptions on the degree of elasticity (i.e. elastic or 
inelastic, never completely inelastic) of the collisions. The 
results on oblique impact physics, reviewed by |Pierazzo fc| 
Melosh (20001, show that oblique impactors ricochet retain- 



ing a significant fraction of the impact speed: if such frac- 
tion is higher than the escape velocity from the target, the 



impactor would bounce and escape into space (Pierazzo & 
Melosh|2000 i. Since the lower limit of the impact velocity we 



found in our calculations is 1km/ s, an order of magnitude 
higher than Phoebe's escape velocity, it is highly probable 
that the impactor's ricochet would allow it to escape from 
the gravitational attraction of the irregular satellite. 
Another important feature of oblique impacts is their lower 
efficiency than head-on impacts in producing high-pressure 



shock waves (Pierazzo & Melosh 20001, making them less 



effective in disrupting the target. Taking into account the 
range of masses and impact velocities which could produce 
Phoebe's capture, this feature of oblique impacts could ex- 
plain Phoebe's surviving to the capture event. 
These results are not meant to be considered as proofs of the 
collisional capture of Phoebe, since the system of equations 
has no unique solution due to the free parameters present 
(i.e. the physical characteristics of the impactor and the 
amount of momentum carried away), but they are never- 
theless the first indication of the viability of the collisional 
capture scenario. Moreover, we found an indirect confirma- 
tion of our results in the analysis of Jason crater performed 



by Gicsc ct al. (2006), who independently observed that the 



morphology of the crater seemed to suggest it originated fol- 
lowing an oblique impact. 

Before proceeding to the next and final section, we need to 
note that, following Gault & Wcdckind] ( |1978[ ), oblique im- 
pacts in the range of values covered by the impact angle in 
both the elastic and inelastic approximations are expected 
to take place with a very low frequency (n w 0.12% in the 
elastic approximation and n « 0.03% in the inelastic one). 
However, the results on Phoebe's collisional efficiency in re- 
moving near-by bodies we presented in Paper I could assist 
us in the interpretation of this issue of the scenario: we will 
discuss the subject in the next section. 



5 CONCLUSIONS 

The aim of this work was to spread a new light on the origin 
of the system of irregular satellites of Saturn and to gain a 
deeper insight on the viability and the consequences of the 
collisional capture process. 

To achieve these goals, we developed the MSSCC software to 
study the inverse capture problem and to search for the pos- 
sible heliocentric pre-capture orbits. Through the use of the 
code, we estimated the specific impulse needed to capture 
the four irregular satellites we concentrated on: Albiorix and 
Siarnaq for the prograde group and Mundilfari and Phoebe 
for the retrograde one. The minimum values of the neces- 
sary impulse we found ranged between 450 — 500 m/s with 
the sole exception of Phoebe which, due to its inner orbit, 
required a minimum value of about 650 m/s. We analysed 
the structure of the solutions in the main orbital elements 
space and we found that those of prograde and retrograde 
satellites were characterised by different features. 
We compared the candidate primordial orbits we computed 
with those of the present day populations of minor bodies 
of the outer Solar System and with the synthetic family of 
the Saturn Crossers found in the original simulations of the 
Nice model. We found that both comets and Centaurs could 
supply candidates to the capture. There was also a good 
match between the synthetic Saturn Crossers and our so- 
lutions, thus supplying the scenario delineated by the Nice 
model with a natural way to populate the system of the 
giant planets with new satellites. We compared the obser- 
vational data available on the irregular satellites with those 
on the minor bodies matching the dynamical range of the 
solutions to the inverse capture problem, and we found that, 
within the errors, there is enough match to strengthen our 
dynamical findings. 

Finally, we used the data supplied by Cassini to test the vi- 
ability of the collisional capture process in non-catastrophic 
impact scenarios and to constrain the impact configuration 
and the impactor mass by comparing our results with the 
information available on Jason crater. We found that highly 
oblique impacts could have produced Jason crater and sup- 
plied the right order of magnitude of specific impulse, the 
requirement on obliquity being imposed by the high specific 
impulse needed to capture Phoebe. The range of values of 
the impact angle associated with our calculations is char- 
acterised by a low probability, making a "single collision" 
scenario to appear implausible. However, in interpreting the 
results exposed we must take into account the following 
fact. Phoebe's surface indicates that the satellite had an 
extremely intense collisional history, with at least 7 craters 
bigger than 50 km and 2 (Jason and Eurytus) of the order 
of 100 km. As a consequence, the dynamical evolution of the 
satellite likely suffered more than one major change in the 
past and the formation of Eurytus crater could be linked 
to Phoebe's capture as well. Since we lack detailed infor- 
mation on Eurytus' morphology, we cannot refine further 
the model. However, a "two collisions" scenario would relax 
the constraints on the impact angle, thus making the whole 
scenario more plausible. Moreover, the results we reported 
in Paper I indicate that Phoebe's efficiency in collisionally 
removing bodies on retrograde orbits is significant only in 
a limited region of phase space (i.e. the one that would be 
occupied by possible collisional fragments). Since we do not 
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observe retrograde satellites in Phoebe's gap (i.e. the region 
comprised between 11.22 x 10 6 km and 14.97 x 10 6 km from 
Saturn, see Paper I) and Phoebe's sweeping effect cannot be 
invoked to explain this fact, the efficiency of the capture pro- 
cess should be low enough to explain the absence of other 
retrograde satellites. The low frequency associated to the 
needed impact parameters in the "single collision" scenario 
could help explaining why Phoebe is the only retrograde 
satellite residing in this orbital region. 
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